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Abstract 

We study discrete vortices in the anti-continuum limit of the discrete two-dimensional non- 
linear Schrodinger (NLS) equations. The discrete vortices in the anti-continuum limit represent 
a finite set of excited nodes on a closed discrete contour with a non-zero topological charge. 
Using the Lyapunov-Schmidt reductions, we find sufficient conditions for continuation and ter- 
mination of the discrete vortices for a small coupling constant in the discrete NLS lattice. An 
example of a closed discrete contour is considered that includes the vortex cell (also known as 
the off-site vortex). We classify the symmetric and asymmetric discrete vortices that bifurcate 
from the anti-continuum limit. We predict analytically and confirm numerically the number of 
unstable eigenvalues associated with various families of such discrete vortices. 

1 Introduction 

Following the first paper of this series 0, we address discrete systems and differential-difference 
equations, which have become topics of increasing physical and mathematical importance. The 
variety of physical applications where such models are relevant, and their significant differences 
from the mathematical theory of partial differential equations, contribute to recent interest in these 
topics. The applicability of such models extends to areas as diverse as nonlinear optics, atomic and 
soft condensed-matter physics, as well as biophysics: specific details and references can be found in 
our first paper P as well as in reviews |2 El E El El El • 

The second paper is devoted to existence and stability of coherent structures in two-dimensional 
lattices, which include both discrete solitons OCI and discrete vortices B, . These two-dimensional 
coherent structures have emerged recently in studies of photorefractive crystals in nonlinear optics 
IHl El and droplets of optical lattices in Bose-Einstein condensates j^fl E] . A significant boost to 
this subject was given by the experimental realization of two-dimensional photonic crystal lattices 
with periodic potentials based on the ideas of [^l. As a result, discrete solitons were observed in 
I14| . while more complex structures such as dipoles, soliton trains and vector solitons were 
observed in |T31[TCl[T7] . 



1 



Most recently, observations of discrete vortices were reported by two independent groups ^1 IT!?] 
where the fundamental vortices with topological charge one were experimentally created and de- 
tected in photorefractive crystals. Two main examples of charge-one discrete vortices include a 
vortex cross (an on-site centered vortex) and a vortex cell (an off-site centered vortex). These 
structures were also recently predicted in a continuous two-dimensional model with the periodic 
potential j2()| . 

These discoveries have stimulated further theoretical work and numerical computations. Thus, 
while in [H], discrete vortices of charge two were shown to be unstable, recently in [23 discrete 
vortices of charge three were found to be stable. Based on the theoretical predictions of |21j . 
further experiments on localized structures were undertaken to unveil other interesting structures, 
such as the discrete soliton necklace which is more globally stable compared to the charge three 
vortex 121] . The discrete vortices have been also extended to three-dimensional discrete models [33] . 
Furthermore, asymmetric vortices have been recently predicted in the two-dimensional lattices in 

El 

The above activity clearly signals the importance and experimental relevance of discrete solitons 
and vortices in two-dimensional discrete lattices. However, most of the above mentioned works are 
predominantly of experimental or numerical nature, while the mathematical theory of existence 
and stability of discrete localized structures has not been developed to a similar extent. The aim 
of the present paper is to develop a categorization of discrete solitons and vortices in the discrete 
two-dimensional nonlinear Schrodinger (NLS) equations. We start from a well-understood limit 
(the so-called anti-continuum case of zero coupling between the lattice nodes) and examining per- 
sistence of the limiting solutions for small coupling by means of the Lyapunov-Schmidt theory. This 
method allows us to discuss persistence and stability of the localized structures by analyzing finite- 
dimensional linear eigenvalue problems. The theoretical predictions agree well with full numerical 
computations of the discrete two-dimensional NLS equation. 

Our main results are summarized for the simplest localized structures in Table 1. These results 
corroborate and extend the previously reported experimental and numerical findings. We quantify 
the stability of the charge-one vortex in accordance to |H] E] El jUI] , the instability of the charge- 
two vortex in accordance to |S1|2] and the stability of the charge-three vortex in accordance to [2]. 
We further demonstrate the instability of all asymmetric vortices proposed in [23]. Furthermore, 
our results can be used to extract the spectral stability of the dipole mode considered in ^5] and 
of the soliton necklace of [22] ■ 
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Table 1: The numbers of small unstable, stable and zero eigenvalues, associated to vortices of the 
discrete two-dimensional NLS equation with small coupling constant. 



The paper is structured as follows. Abstract results on existence of discrete solitons and vortices 
are derived in Section 2. Persistence of localized modes for a particular square discrete contour is 
considered in Section 3. Stability of the persistent solutions is addressed in Section 4. Analytical 
results are compared to numerical computations in Section 5. Section 6 concludes the paper with 
a summary of our results and discussions of interesting directions for future study. 



2 Existence of discrete vortices 

We consider the discrete nonlinear Schrodinger (NLS) equation in two space dimensions 

where Un,m{t) ■ IR+ C, {n,m) S 1?, and e > is the inverse squared step size of the discrete 
two-dimensional NLS lattice. Time-periodic localized modes of the discrete NLS equation 1)2. l|l take 
the form: 

=0n,me*('^-4^)*+*^°, 0„,™eC, (n,TO)eZ2, (2.2) 

where G and /i G M arc parameters. Since localized modes in the focusing NLS lattice (12.1(1 with 
e > may exist only for /i > 4e 4 and the parameter /i is scaled out by the scaling transformation, 

4>n,m ^ \/T^4'n,rn, £ = M^, (2.3) 

the parameter > will henceforth be set to = 1. In this case, the complex- valued 0„,m solve 
the nonlinear difference equations on {n,m) E 1?: 

(1 - \<i>n,m}\')<\>n,m = £ [<\>n+\,m + 0ri-l,m + <i>n,m+\ + ^n,rn-\) ■ (2.4) 

As e = 0, the localized modes of the difference equations 1)2.4(1 are given by the limiting solution: 

,(0) _ f e^'""", {n.m) e 5, 
0, (n,m)eZV, 

where S" is a finite set of nodes on the lattice (n, to) G 1? and Qn.m are parameters for (n, to) G S. 
Since 0o is arbitrary in the ansatz 1(2.2(1 . we can set On„,mo = for a particular node (no,TOo) G S. 
Using this convention, we define two special types of localized modes, called discrete solitons and 
vortices. 

Definition 2.1 The localized solution of the difference equations \2.4\j with e > 0, that has all real- 
valued amplitudes 4>n,m> (n^m) G and satisfies the limit 1^2. 5\} with all 9n^m = {0, tt}, (n,m) G S, 
is called a discrete soliton. 

Definition 2.2 Let S be a simply- connected discrete contour on the plane {n,m) G Z^. The local- 
ized solution of the difference equations {2.4-^ with e > 0, that has complex-valued 4>n,m, {n,m) G Z^ 
and satisfies the limit 1^2. 5\) with 9n,m G [0, 27r], {n,m) G S, is called a discrete vortex. 
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Definition 2.3 Let S be a simply- connected discrete contour on the plane {n,m) G 1? , such that 
each node (n, m) G S has exactly two adjacent nodes in vertical or horizontal directions along S . If 
the phase difference between two adjacent 9n,m for {n, m) G S is constant in S, the discrete vortex 
is called symmetric. Otherwise, it is called asymmetric. The total number of 2t: phase shifts across 
the closed contour S is called the topological charge of the discrete vortex. 

In particular, we consider the ordered simply-connected discrete contour S = Sm: 

Sm = {(1, 1), (2, 1), (M + 1, 1), (M + 1, 2), {M + 1, M + 1), 

(M, M + 1), (1, M + 1), (1, M), (1, 2)} , (2.6) 

where dim(5'Af) = 4M. According to Definition 12.31 the contour Sm for a fixed M could support 
symmetric and asymmetric vortices with some charge L. An example of the simplest vortices for 
M = 1 is the symmetric charge-one vortex cell (6*1,1 = 0, 6*2,1 ~ f , &2.2 ~ tt, ^1,2 = ^) [HI HO] and 
an asymmetric charge-one vortex (6*1,1 = 0, 62,1 — 0, 62,2 — tt, 6*1,2 = ti" + 6*), where 6* ^ {0, f ,7r} 

m 

It follows from the general method [3 ESI that the discrete solitons of the two-dimensional NLS 
lattice (|2.4|l (see Definition 12.1(1 can be continued to the domain < e < eq for some eo > 0. 
It is more complicated to find a configuration of 9n,m for (n, m) S S that allows us to continue 
the discrete vortices (see Definition 12.211 for e > 0. The continuation of the discrete solitons and 
vortices is based on the Implicit Function Theorem and Lyapunov-Schmidt Reduction Theorem 
[261 177] . Abstract results on existence of such continuations are formulated and proved below, after 
the introduction of some relevant notations. 

Let 0(0) be a small neighborhood of e = 0, such that 0(0) = (— eo, ^o) for some eo > 0. Let 
= dim(S') and T be the torus on [0, 2tt]^ , such that 6n.m for {n, m) e 5 form a vector 9 G T. Let 
J7 = L^(Z^, C) be the Hilbert space of square-summable complex-valued sequences {4>n,m}{n,m)£Z^, 
equipped with the inner product and the norm: 

(U,v)n= ^ Un,mVn,7n, I|u||i2= ^ \Un,m\'^<00. (2.7) 

Let u denote an infinite-dimensional vector in Q, that consists of components Un.m for all {n, m) G Z^. 

Proposition 2.4 There exists a unique (discrete soliton) solution of the difference equations \2.4^ 
in the domain e S 0(0) that satisfies (i) S K, {n,m) G J? and (ii) limg_+o = (t)n}m, where 
(l)n}m is given by 1^2. 5\) with O^.m — {0,7r}, {n,m) G S. The solution 4>{e) is analytic in e G 0(0). 

Proof. Assume that 4)n,m G M for all {n,m) G 1?. The difference equations (|2.4|l are rewritten as 
zeros of the nonlinear vector-valued function: 

fn,m{4>, e) = (1 - (t)n^,n)4>n,m ~ e {4>n+l,m + 4>n-l,m + 4>n,m+l + 0n,m-l) = 0. (2.8) 

The mapping f : x 0(0) i-^ is on G and has a bounded continuous Frechet derivative, 
given by: 

Cn,ra = (l - 30^, - e (s + i,o + 5-1,0 + S0, + 1 + SO,-l) ■ (2.9) 

where Sn',m' is the shift operator, such that Sn',m'Un.m = Un+n'.rn+m'- It is obvious that 

f(0("\O) = O, ker(/:(°)) = 0, (2.10) 
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where 4>^^^ is the discrete sohton of Definition 12 . II and ^C*^"^ is the operator C computed at = 
and e = 0. ft follows from (|?l^ - (|?lIHl that : 17 i-> has a bounded inverse. By the Implicit 
Function Theorem 1371, Appendix f], there exists a local mapping cj) : 0(0) — > fl, such that ^(e) 
is continuous in e e O{0) and 4>^^'' — 0(0). Moreover, since f(0, e) is analytic in e e 0(0), then 
0(e) is analytic in e e 0(0) ^ Chapter 2.2]. ■ 



Remark 2.5 Proposition \2.Ji\ does not exclude a possibility of continuation of the limiting solution 
1^2. 5\} with 9n,m — {0,7r} for all (n,m) & S to the complex-valued solution <p{e) in e £ 0(0). 

Proposition 2.6 There exists a vector-valued function g : T x 0(0) ^-^ M^, such that the limiting 
solution \2.5]) is continued to the domain e € 0(0) if and only if 6 € T is a root of g{6,e) — in 
e e 0(0). Moreover, the function g(0, e) is analytic in e € 0(0) and g{6, 0) = for any 6 (z T. 

Proof. When (f)n.m G C for some {n, m) £ 1? , the difference equations H2.4|l are complemented by 
the complex conjugate equations in the abstract form: 

f(0,0,e) = O, f(0,0,e) = O. (2.ff) 

Taking the Frechet derivative of f(0, 0, e) with respect to and 0, we compute the linearization 
operator H for the difference equations (|2.4|) : 

Hn,™ = f ^"^K"'™'' ^_~tr ,2 V^(^+M+5-i,o+5o,+i+=so,-i)i : : ). (2.12) 




Let 'H^°'> ^ 7i:(0'°\ 0). ft is clear that H^°'> : nxn^nxftisn self-adjoint Fredholm operator of 
index zero with dim ker(7i'^*'^) — N. Moreover, eigenvectors of ker(7i^°)) re-normalize the parame- 
ters On,m for [n,m) G 5 in the limiting solution p. 5(1 . By the Lyapunov Reduction Theorem |27l 
Chapter 7.1], there exists a decomposition 17 = ker(?i'^'^') © w, such that g(0, e) is defined in terms 
of the projections to ker(7i*^°''). Let {e„^m}(n,m)GS be a set of N linearly independent eigenvectors 
in the kernel of Ti'"^. ft follows from the representation, 

= - ^_2l„,„ ^''i"'" ) ' (^^'™) e (2.13) 

that each eigenvector e„_m in the set {&n,m}{n,m)^s has the only non-zero element (e*^" ™, — e"*^"."")"^ 
at the (n,m)-th position of u G 12. By projections of the nonlinear equations (|2.11l) to ker(7i'^''^), 
we derive an imphcit representation for the functions g(0, e): 

2ign,m{0, e) = (1 - l^n.mP) (e~*^"'"0„,,„ - e'^"-'"0„,™) 

n+l,m ) , (2.14) 

for (n,m) £ 5, where the factor (2i) is introduced for convenient notations. Let (j>n,m = e'^" '"'itn,m 
for {n,m) € S and 0„^m = for {n,m) £ 1?\S. Since eigenvectors of ker(H^°-') are excluded 

from the solution in w C 17, we have Un,m € R for ("■, ™) G 5*, such that 

- 2ign,miO,e) = ee""^"'" (0„+l,„i + 0n-l,m + 0n,m+l + 0n,m-l) 

- ee*^"'" (0„+l,„i + 0„-l,m + 0n,m+l + 0„,m-l) (2-15) 
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and g(0, 0) = for any €T. Since f(</>, e) is analytic in e G 0(0), then g(0,e) is analytic in 
e e 0(0) 27, Appendix 3]. ■ 

Corollary 2.7 T/ie Junction g(0, e) can 6e expanded into convergent Taylor series in 0(0); 

oc 

g(0,6)=^e'^gW(0), g«(0) = -a,^g(0,O). (2.16) 

fc=i 

IJ the root 6{e) ofg{9,e) ~ is analytic m e G 0(0), then the solution 0(e) is analytic m e G 0(0), 
such that 

oc 

0(e) = 0("'+5]eV\ (2.17) 



fc=i 



where (fr ' is given by fij.,5|) . 

Lemma 2.8 Let 0(e) 6e a root of g{6,e) ~ in e E 0(0). An arbitrary shift 6{e) + 6'oPo, where 
6*0 G K and po = (1,1,..., 1)^ , gives a one-parameter family of roots of g(0, e) = for the same e. 

Proof. The statement follows from the symmetry of the difference equations (|2.4|l with respect 
to gauge transformation Tf, Chapter 7.3]. ■ 

Proposition 2.9 Let 9^ be the root of ^^^\0) ~ and Aii be the Jacobian matrix o/g^^^(0) 
at 9 ^ 9^. If the matrix M.i has a simple zero eigenvalue, there exists a unique (modulo gauge 
transformation) analytic continuation of the limiting solution \2.5]) to the domain e G 0(0). 

Proof. By Lemma 12.81 the matrix Ai i has always a non-empty kernel with the eigenvector 
Po = (1, 1, 1) due to gauge transformation. Let be the constrained subspace of C^: 

Xo = {ugC": (po,u) = 0}. (2.18) 

If the matrix A^i is non-singular in the subspace Xq, then there exists a unique (modulo the shift) 
analytic continuation of the root 0* in e G 0(0) by the Implicit Function Theorem, applied to the 
nonlinear equation g(0, e) = [371 Appendix 1]. ■ 

Proposition 2.10 Let 9^ 6e a (1 -f- d)-parameter solution of g^^^ (9) = and A4i have a zero 
eigenvalue of multiplicity (1 -I- d), where 1 < d < N — 1. The limiting solution 1^2. 5\) can be 
continued in the domain e G 0(0) only if g^^^{9^,) is orthogonal to ker(A^i). 

Proof. Let po and {pzjf^]^ be eigenvectors of ker(A^i). We define the constrained subspace of 
Xo: 

Xd^{ueXo: (pi,u) = 0, ; = l,...,d}. (2.19) 

If g'^^^(0*) ^ Xd, the Lyapunov-Schmidt Reduction Theorem in finite dimensions |27l Chapter 1.3] 
shows that the solution 9^ can not be continued in e G 0(0). ■ 

Prot)osition l2.6l gives an abstract formulation of the continuation problem for the limiting solution 
(|2.5|l for e 7^ 0. ProDOsition l2 . 91 gives a sufhcient condition of existence and uniqueness (up to gauge 
invariancc) of such continuations. Proposition 12. 1(11 gives a sufhcient condition for termination of 
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multi-parameter solutions. Particular applications of Propositions 12 . til 12 . 91 and 12 . 1 Ul are limited by 
the complexity of the set S in the limiting solution H2.5|l . since computations of the vector- valued 
function g*^^)(0), g(^'(0), and the Jacobian matrix Aii could be technically involved. We apply 
the abstract results of Propositions l2.6l ITHI and l^.lUl to the simply-connected discrete contour Sm, 
defined in 



3 Persistence of discrete vortices 



We consider discrete solitons and vortices on the contour Sm defined by H2.6|l . Let the set 9j 
correspond to the ordered contour Sm, starting at 9i = 9i^i, 02 = 02, i and ending at 9n = ^1,2, 
where N = AM . In what follows, we use the periodic boundary conditions for 9j on the circle from 
j 1 to j = N, such that 6*0 ~ 9^, 9i = ^^w+i, and so on. 

The discrete vortex has the charge L if the phase difference changes on 27rL along the discrete 
contour By gauge transformation, we can always set 9i = for convenience. We will also 

choose 92—9 with < < tt for convenience, which corresponds to discrete vortices with i > 0. 



3.1 Solutions of the first-order reductions 

Substituting the limiting solution (j)n.m, in the bifurcation function H2.15(l . we find that g'-^\e) m 
the Taylor series H2.16|l is non-zero for the contour 5a/ and it takes the form: 

gf\e) = sm{9j - 9^+,) + sin(0, - 0,_i), 1 < j < A^. (3.1) 

The bifurcation equations g'^^^(0) = are rewritten as a system of N nonlinear equations for N 
parameters 9i,92,---,9n as follows: 

sin(02 - Oi) = sin(03 - 02) - ... = sm{9N - 9n-i) = sm{9i - 9n). (3.2) 

We classify all solutions of the bifurcation equations (|3.2|l and give explicit examples for M = 1 
and M = 2. 

Proposition 3.1 Let aj = cos{9j+i — 9j) for I < j < N , such that 9i = 0, 92 = 9, and ^at+i — 2ttL, 
where N ~ AM , < 9 < n and L is the vortex charge. All solutions of the bifurcation equations 
ij^*.^) reduce to the four families: 

(i) discrete solitons with 9 = {0, tt} and 

^^,={0,7r}, 1<3<N, (3.3) 



< / < iV. 



such that the set {ojjjLi includes I coefficients aj — 1 and N — I coefficients aj = —I, where 



(a) symmetric vortices of charge L with 9 — where 1 < L < 2M — 1, and 
such that all N coefficients are the same: aj = a — cos (^^) . 
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(Hi) one-parameter families of asymmetric vortices of charge L = M with < 9 < n and 

Oj+i - = I I mod(2^), 2<j<N, (3.5) 

such that the set {ajJ^L^ includes 2M coefficients aj — cos6 and 2M coefficients aj — — cos0. 
(iv) zero-parameter asymmetric vortices of charge L ^ M and 



n-2M 



1 < n < TV - 1, n^2M, (3.6) 



such that the set {ajJ^L^ includes n coefficients aj — cos0* and N — n coefficients aj = — cos0* 
and the family (iv) does not reduce to any of the families (i)-(iii). 

Proof. All solutions of the bifurcation equations H3.2|l are given by the binary choice (|3.5|l in the 
two roots of the sine-function on 6* € [0, 27r], where the first choice gives Oj = cos 6* and the second 
choice gives Oj — — cos 0. Let us assume that there are totally n first choices and N — n second 
choices, where 1 < n < A^. Then, we have 

On+i ^ne + {N - n){'K -9) = {2n - N)9 + {N - n)n = 2nL, 

where L is the integer charge of the discrete vortex. There are only two solutions of the above 
equation. When 9 is arbitrary parameter, we have n — ^ — 2M and L — M, which gives the 
one-parameter family (iii). When 9 — 9^, is fixed, we have 

TT f n + 2L - AM 
n - 2M 

When n — N — 2L, we have the family (i) with N — 2L phases 9j — and 2L phases 9j — ir. 
Since discrete solitons do not have topological charge, the parameter L could be half-integer: L — 
{N - 0/2, where < I < N. When n = AM, we have the family (ii) for any 1 < L < 2M - 1. 
Other choices of n, which are irreducible to the families (i)-(iii), produce the family (iv). ■ 

Remark 3.2 The one-parameter family (iii) connects special solutions of the families (i) and (ii). 
When 9 = and 9 = the family (iii) reduces to the family (i) with I = 2M . When 9 = ^, the 
family (iii) reduces to the family (ii) with L = M. We shall call the corresponding solutions of 
family (i) as the super- symmetric soliton and of family (ii) as the super- symmetric vortex. 

Remark 3.3 There exist Ni = 2^^-'^ solutions of family (i), N2 = 2M ~ 1 solutions of family (ii), 
and solutions of family (iii), where 

2M-1 
fc=0 ^ ' 

The number A^4 of solutions of family (iv) can not be computed in general. We consider such 
solutions only in the explicit examples of M = 1 and M = 2. 

Example M = 1 and = 4: There are A^i = 8 solutions of family (i), A^2 = 1 solution of family 
(ii), = 3 solutions of family (iii), and no solutions of family (iv). The three one-parameter 
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asymmetric vortices are given explicitly by 

(a) 01 = 0, 02 

(b) 01 ^ 0, 02 

(c) 01 - 0, 02 



^0, 03 = TT, 04 = 7r + (3.8) 
= 0, 03 = 20, 04 TT + (3.9) 

= 0, 03 = TT, 04 = 27r - 0. (3.10) 



Example M = 2 and = 8: There are A'^i = 128 solutions of family (i), A^2 = 3 solutions 
of family (ii), A^3 — 35 solutions of family (iii), and A'4 — 14 solutions of family (iv). The three 
symmetric vortices have topological charge L = 1 {0 — j), L ^ 2 {0 ^ and L — 3 {0 — ■^). The 
one-parameter asymmetric vortices include 35 combinations of 4 upper choices and 4 lower choices 
in (|3.5|) . starting with the following three solutions: 

(a) 01 = 0, 02 = 0, 03 = 20, 04 = 30, 05 = 40, 06 = TT + 30, 07 = 27r + 20, 08 = 37r + 0, 

(b) 01 = 0, 02 = 0, 03 = 20, 04 = 30, 05 = TT + 20, 06 = TT + 30, 07 = 27r + 20, 08 = 37r + 0, 

(c) 01 = 0, 02 = 0, 03 = 20, 04 = 30, 05 = TT + 20, 06 = 27r + 0, 07 = 27r + 20, 08 = 37r + 0, 

and so on. The zero-parameter asymmetric vortices include 7 combinations of vortices with L = 1 
for seven phase differences ^ and one phase difference ^ and 7 combinations of vortices with L = 3 
for one phase difference | and seven phase differences 

3.2 Continuation of solutions of the first-order reductions 

We compute the Jacobian matrix Aii from the bifurcation function g'^^^(0), given in (|3.1|l : 

{cos(0j+i - 0j) + cos(0j_i - 0j), i = j, 

-cos(0j-0,), i^j±l (3.11) 

0, K-j|>2 

subject to the periodic boundary conditions. The matrix Mi is defined by the coefficients aj = 
cos(0j-|_i — 0j) for 1 < j < A^. It has the same structure as that in the perturbation theory of 
continuous multi- pulse solitons in coupled NLS equations |2H]. Three technical results establish 
location of eigenvalues of the matrix Aii. 

Lemma 3.4 Let uq, zq, and po be the numbers of negative, zero and positive terms of Oj = 
cos(0j-|_i — 0j), 1 < j < N , such that uq + zq + po = N. Let n{A4i), z{A4i), and p{A4i) be 
the numbers of negative, zero and positive eigenvalues of the matrix Mi, defined by 15*.^ Assume 
that zq = and denote: 

^^=tl[^.= (f[A (E^)- (3-12) 

1=1 \i=l J \i=\ V 

If Ai ^ 0, then, z{Mi) = 1, and either n{Mi) = uq — 1, p{Mi) = po or n{Mi) = uq, p{Mi) — 
Pq — 1. Moreover, n(Mi) is even if Ai > and is odd if Ai < 0. If Ai = 0, then z(A4i) > 2. 

Proof. The first statement follows from Appendix A of [211. Let the determinant equation be 
D{\) = det(A^i — XI) = 0. By induction arguments in ISHIEHIj it can be found that D{0) — and 
D'{0) — —NAi. On the other hand, -D'(O) = — A1A2 • • • Aat-i, where Xn = (which exists always 
with the eigenvector po = (1, 1, 1)^, see Proposition 12. 9|l . Then, it is clear that (— l)"(-^i) = 
sign(Ai). When Ai = 0, at least one more eigenvalue is zero, such that z{Aii) > 2. m 
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Lemma 3.5 Let all coefficients Uj = cos(0j+i — 9j), 1 < j < N be the same: aj = a. Eigenvalues 
of the matrix A4i are computed explicitly as follows: 

A„=4asin^^, l<n<N. (3.13) 

Proof. When aj = a, 1 < j < N, the eigenvalue problem for the matrix A4i takes the form of 
the linear difference equation with constant coefficients: 

a {2xj — Xj^i — Xj-i) = Xxj, xq = xn, xi — xn+i, (3-14) 

The discrete Fourier mode Xj = exp (? ^^" ) for 1 < j,n < N results in the solution (|3.13l) . ■ 

Lemma 3.6 Let all coefficients aj = cos{9j+i — 9j), 1 < j < N alternate the sign: aj ~ (— l)-'a, 
where N = AM . Eigenvalues of the matrix Mi are computed explicitly as follows: 

TTTl 

A„ = -A„+2Af = 2asin— , 1 < n < 2Af, (3.15) 

such that n{Mi) — 2M — 1, z{Mi) — 2, and p{Mi) — 2M — 1. These numbers do not change if 
the set {aj}jLi is obtained from the sign- alternating set by permutations. 

Proof. When aj = (— l)-'a, I < j < AM, the eigenvalue problem for the matrix Aii takes the 
form of a coupled system of linear difference equation with constant coefficients: 

O'iyj ~ Vj-i) ^ ^Xj^ a{xj - Xj+i) ^ Xyj, 1 < J < 2M, (3.16) 

subject to the periodic boundary conditions: xi — a;2_A/+i and j/o = 2/2Af ■ The discrete Fourier mode 
Xj = Xq exp {i^§^) and yj = yo exp {i^§^) for 1 < j,n < 2M results in the solution (|3.15|) . In this 
case, we have n{Mi) = 2M - 1, z{Mi) = 2, and p{Mi) = 2M - 1, such that ^(O) = D'{0) = in 
the determinant equation D{X) = det(A^i — A/). In order to prove that z{Mi) — 2 remains invariant 



with respect to permutations of the sign-alternating set {aj}jLi, we find from Mathematica that 

D"{0) ^ (f[ a) (a^ + + f E E ^ -] ] ■ (3-17) 

where < < Pn are numerical coefficients. Let A* denote the sign-alternating set {ajlj^i, 
such that Oj ~ (— l)-'a, and A denote a set obtained from by permutations. It is clear that 

and 

/N-l N \ /JV-1 N 



\ 1=1 1=1+1 'A, ^ /=J-I-1 / ^ 

Therefore, the expression in brackets in H3.17|l can be estimated as follows: 

\ 1=1 /A \ '=1 i=2+l / A 

where the last inequality follows from the fact that £'"(0) < for A*. Therefore, z{Aii) = 2 for A. 
Combining it with estimates from Appendix A in |2H|, we have n{A4i) = p{Aii) = 2M — 1. ■ 
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Using Leinmas l3.4ll3.5l and l3.5l we classify the continuation of solutions of the first-order reductions, 
which are described in the families (i)-(iv) of Proposition 13. II 

For family (i), excluding the case of super-symmetric solitons (see Remark \'6.'2\i . the numbers of 
positive and negative signs of Uj are different, such that the conditions 2:0 = and Ai ^ are 
satisfied in Lemma [3.41 and hence z{A4i) = 1. By Proposition 12.91 the family (i) has a unique 
continuation to discrete solitons (see Definition 12. 1|) . Continuations described in Remark 12.51 are 
only possible for super-symmetric solitons. 

For family (ii), all coefficients aj are the same: aj = a = cos 

(jm)^ < j < By Lemma ESI 
there is always a zero eigenvalue (Xn =0), while remaining {N — 1) eigenvalues are all positive for 
a > (when 1 < L < M - 1), negative for a < (when M + 1 < L < 2M - 1), and zero for a = 
(when L = M). By Proposition l2.9l the family (ii) has a unique continuation to symmetric vortices 
with charge L, where 1 < i < 2Af - 1 and M (see Definitions lO and lO) . 

For family (iii), there are 2M coefficients Uj — cos 9 and 2M coefficients Uj = — cos 9, which are 
non-zero for 6* ^ f . By Lemma IX^ we have n{Mi) = 2M - 1, z{Mi) = 2, and piMi) = 2M - 1. 
The additional zero eigenvalue is related to the derivative of the family of the asymmetric discrete 
vortices (|3.5|) with respect to the parameter 9. Therefore, continuations of family (iii) of asymmetric 
vortices, including the particular cases of super-symmetric solitons of family (i) and super-symmetric 
vortices of family (ii), must be considered beyond the first-order reductions. 

For family (iv), since n ^ 2M, the conditions zq — and Ai ^ are satisfied in Lemma f3.4l and 
hence z{A4i) = 1. By Proposition 12.91 the family (iv) has a unique continuation to asymmetric 
vortices for e 7^ 0. 

3.3 Continuation of solutions to the second-order reductions 

Results of the first-order reductions are insufficient to conclude persistence of the asymmetric vor- 
tices of family (iii), including the super-symmetric soliton of family (i) and the super-symmetric 
vortex of family (ii). Therefore, we continue the bifurcation function g(0,e) to the second order of 
e in the Taylor series 12.16|l . It follows from 12.4(1 that the first-order correction of the Taylor series 
(|2.17|l satisfies the inhomogeneous problem: 

(1 - 2i<Ln0ii - cf>i^mL = €i,m + + <Ui + (3-18) 

We define solution of the inhomogeneous problem 1)3.18(1 in w C 51, such that the homogeneous 
solutions in ker(7i(°)) are removed from the solution (p'^^K This is equivalent to the constraint: 
0n,m — u„.„ie*^" '" , Un.m S R for all {n,m) £ Sm- We develop computations for three distinct 
cases: M — 1, M = 2 and M > 3. This separation is due to the special structure of the discrete 
contours Sm- 

Case M = 1: The inhomogeneous problem 1(3.18(1 has a unique solution t/)^^-' G C 51: 

= [cos((?,_i - 6,) + cos(0,+i - 6,)] e*«^ (3.19) 
where the index j enumerates the node (n, m) on the contour Sm, 

</'il = e^'^ (3.20) 
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where the node (n, m) is adjacent to the j-th node on the contour Sm, while 



is empty for all 



into the bifurcation function 



remaining nodes. By substituting the first-order correction term (l>n)rn 
H2.15|l . we find the correction term g'^^'>{9) in the Taylor series (|2.16|l : 

ij) [cos{9j - 9j-i) + cos{9j-2 - Oj-i)] , 1 < i < ^. 
We compute the vector g2 = g'--^\0) at the asymmetric vortex solutions (|3.8|) - (|3.1U|) : 



- sm(t/j+i 
isin(6lj_i 



(3.21) 



(a) g2 








(b) g2 



/ 2 \ 



-2 

V J 



sm f cos ( 



(c) g2 



f \ 

-2 



V 2 y 



sm f cos ( 



The kernel of A^i is two-dimensional with the eigenvectors po and pi. The second eigenvector pi 
is related to derivatives of the solutions H3.8|l - H3.10|l in 9: 



(a) Pi = 



1 


V 1 y 



(b) Pi = 



1 
2 

V 1 y 



(c) Pi = 



f \ 

1 



V -1 y 



The Fredholm alternative (pi,g2) = is satisfied for the solution (a) but fails for the solutions (b) 
and (c), unless 9 = {0, The latter cases are included in the definitions of super-symmetric 

discrete solitons and vortices (see Remark l3 .21) . By Proposition 12. 101 the solutions (b) and (c) can 
not be continued in e 7^ 0, while the solution (a) can be continued up to the second-order reductions. 

Case AI = 2: The solution (p^^^^ e C of the inhomogeneous problem H3.18|) is given by (|3.19|) 
and (|3.2(J|) . except for the center node (2,2), where 



4!^ 



f2 +g»y4 _,_g^W6 _^g««8_ (3 22) 

The correction term g'^^) (g) given by H3.21|l but the even entries are modified as follows: 



sm f — I 



sm{9j-9j+4), J = 2,4,6,8. (3.23) 



The vector g2 = g'^-* (9) can be computed for each of 35 one-parameter asymmetric vortex solutions, 
starting with the first three solutions: 



(a) g2 = 



/ 2 \ 

1 


-1 
-2 
-1 


V 1 y 



sinf cost 



(b) g2 = 



/ 2 \ 

1 

-1 
-1 


-1 
-1 

V 1 y 



sm cos I 



(C) g2 = 



/ 2 \ 

1 

-1 
-2 

1 

-1 

V y 



sm W cos f 
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The second eigenvector pi of the kernel of A^i is related to derivatives of the family in 9, e.g. 



(a) Pi = 



/ \ 
1 
2 
3 

4 
3 
2 

V 1 / 



(b) Pi 



/ \ 

1 
2 
3 
2 
3 
2 

V 1 / 



(c) Pi 



/ \ 
1 
2 
3 
2 
1 
2 

V 1 / 



The Fredholm alternative condition (pi,g2) — fails for all solutions of family (iii) but one, 
excluding the special values 9 = {0, f,7r}. The only solution of family (iii), where g2 = 0, is 



characterized by the alternating signs of coefficients aj = cos( 



,) for 1 < J < A^. 



Case M > 3: The solution e ti^ C 51 of the inhomogeneous problem H3.18|) is given by (|3.19fl 
and ()3.2l)fl . except for the four interior corner nodes (2, 2),(Af, 2),(Af, M), and (2, Af), where 



j = 1, M + 1, 2M + 1, 3M + 1. 



(3.24) 



The correction term g^^^(0) is given by (|3.21|) . except for the adjacent entries to the four corner 
nodes on the contour ^Af : (1, 1), (1, M + 1), (M + 1, M + 1), and (M + 1, 1), which are modified 
by: 



■'3 



\0) 



sm 
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J+2), 



j ^2,M + 2,2M + 2, 3M + 2, 
j = M, 2M, 3M, 4M. 



(3.25) 



Again, there is only one solution of family (iii), where g2 = 0, which is characterized by the 
alternating signs of coefhcients aj — cos{9j+i — 9j) for 1 < j < A^. All other solutions of family (iii) 
do not satisfy the Fredholm alternative condition (pi,g2) ~ 0. 

By using results of these computations, we classify continuations of solutions of the super-symmetric 
solitons of family (i) and asymmetric vortices of family (iii). Let 7M2 be the Jacobian matrix 
computed from the bifurcation function g^^)^^)^ given in H3.21|l . H3.23|l and H3.25|) . Since (pi, g2) 7^ 
for 9 7^ {0, f-jTr}, except for the case of sign-alternating set {aj}jLi with aj = (— l)^a, it follows 
from regular perturbation theory that (pi,A^2Pi) 7^ 0. Therefore, the second zero eigenvalue of 
A4i bifurcates off zero for the matrix A4i + eM.2- By Proposition |^| (which needs to be modified 
for the Jacobian matrix A^i +eA^2), the super-symmetric solutions of family (i), which are different 
from sign-alternating sets aj = (— l)^a, are uniquely continued to discrete solitons (see Definition 
EH). 

By Proposition 12.101 all asymmetric vortices of family (iii) , except for the sign-alternating set 
aj = cos{9j+i — 9j) — (—1)-'+^ cosO, I < j < N , can not be continued to e 0. The only solution 
which can be continued up to the second-order reductions has the explicit form: 



94j_3 = 27r(i-l), 



Uj-2 



Hj-1 



%j-3 + TT, 



hj-3 + TT + I 



(3.26) 



where 1 < j < M and < 9 < tt. This solution includes two particular cases of super-symmetric 
solitons of family (i) for 9^0 and 9 = tt and super-symmetric vortices of family (ii) for 9 — ^. 
Continuation of the solution 13.26|l must be considered beyond the second-order reductions. 
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3.4 Jacobian matrix of the second-order reductions 



The Jacobian matrix Mi oi the first-order reductions is empty for super-symmetric vortices of family 
(ii) with L = M. In order to study stabihty of super-symmetric vortices, we need to compute the 
Jacobian matrix M.2 from the second-order bifurcation function g^^-'(&), given in (|3.21|l . 13. 2311 . and 
(|3.25|l . These computations are developed separately for three cases M = 1, M = 2, and M > 3. 

Case M = 1: Non-zero elements of A^2 are given by: 



0, 



« = J, 

i=j±2 
\i-j\^0,2 



(3.27) 



or explicitly for = 4: 



M2 = 



( 1 





-1 





\ 





1 





-1 




-1 





1 







V 


-1 





1 





(3.28) 



The matrix Ai2 has four eigenvalues: Ai = A2 = 2 and A3 = A4 = 0. The two eigenvectors for 
the zero eigenvalue are pa = (1, 0, 1, 0)^ and p4 = (0, 1, 0, 1)^. The eigenvector p4 corresponds to 
the derivative of the asymmetric vortex H3.8(l with respect to parameter 9, while the eigenvector 
Po = P3 + P4 corresponds to the shift due to gauge transformation. 

Case M — 2: The Jacobian matrix AI2 is given in H3.27fl except for the even entries which are 
modified as follows: 



-1, «=j, 

iM2)^,J ^ iM2)^,j + { +1, i=j±2 

-1, i=j±A 
The matrix M.2 for = 8 takes the explicit form: 



J = 2, 4, 6, 8. 



(3.29) 



Mo 



1 





1 

2 











1 

2 














1 

2 





-1 





1 

2 


1 

2 





1 





1 

2 














1 

2 











1 

2 





-1 








1 

2 





1 





1 

2 








-1 





1 

2 











1 

2 


1 

2 











1 

2 





1 








1 

2 





-1 





1 

2 









(3.30) 



The eigenvalue problem for M2 decouples into two linear difference equations with constant coef- 
ficients: 



~2yj+2 + yj+i + yj-i 



- 2\xj, 



J 



= 1,2,3,4 
1,2,3,4, 



subject to the periodic boundary conditions for Xj and yj. By the discrete Fourier transform, see the 
proof of Lemma lT^ the first problem has eigenvalues: Ai = 1, A2 — 2, A3 = 1, and A4 = 0, while the 
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second problem has eigenvalues: A5 = 1, Ae = —2, A7 = 1, and Ag = 0. The two eigenvectors for the 
zero eigenvalue are P4 = (1, 0, 1, 0, 1, 0, 1, 0)'^ and ps = (0, 1, 0, 1, 0, 1, 0, 1)"^, where the eigenvector 
Ps corresponds to the derivative of the asymmetric vortex H3.26|l with respect to parameter and 
the eigenvector Po = P4 + Ps corresponds to the shift due to gauge transformation. 

Case AI > 3: The Jacobian matrix A^2 is given in H3.27|l . except for the adjacent entries to the 
four corner nodes on the contours 5*^/: (1, 1), (1, M + 1), (Af + 1, M + 1), and (M + 1, 1), which 
are modified by 



(A^2).j {M2)^^J + { +1, 

fl, 



i = j = 2,M,M + 2, 2M, 2M + 2, 3M, 3M + 2, 4M, 

i = j -2 = M, 2M, 3M, 4:M 

i = j + 2 = 2,M + 2, 2M + 2, 3M + 2 



(3.31) 



The matrix A42 for TV = 12 takes the explicit form: 



M2 



V 



_i 
2 






1 

2 











_1 

2 



1 



_1 

2 














1 

2 





_ 1 

2 













_1 

2 






1 

2 
















2 



1 



_1 

2 



_ 1 

2 










1 

2 







\ 

1 

2 









2 






(3.32) 



The eigenvalue problem for M2 decouples into eigenvalue problems for two 6-by-6 matrix, which are 
related by the Toeplitz transformation. As a result, the spectra of these two matrices are identical 
with the eigenvalues, obtained with the use of MATLAB: 

Ai = A7 = -0.780776, A2 = Ag = -0.5, A3 = Ag = 0, 
A4 = Aio = 0.5, As = All = 1.28078, Ag = A12 = 1.5. 



We confirm that the matrix AI2 has exactly two zero eigenvalues, one of which is related to the 
derivative of the asymmetric vortex H3.26|l in 9 and the other one is related to the shift due to gauge 
transformation. 

Computations of the matrix AI2 for super-symmetric vortices of family (ii) confirm the results of 
the second-order reductions for asymmetric vortices of family (iii). Although all N3 solutions of 
family (iii) reduce to the super-symmetric vortex of family (ii) in the first-order reductions, it is 
the only family H3.26|l that survives in the second-order reductions, such that the super-symmetric 
vortex of family (ii) with L = M and ^ ^ can be deformed and continued up to the second-order 
reductions to the asymmetric vortex H3.26|l . 

All individual results on persistence of localized modes on the discrete contour Sm are summarized 
as follows. 
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Proposition 3.7 Consider the discrete soliton and vortices of the nonlinear equations {2.4^ that 
bifurcate from the limiting solution 4)n]m in \2.5\) on the discrete contour in \2.(^) . There exists 
a unique (modulo gauge transformation) continuation to the domain e € 0{Q) of discrete solitons 
of family (i) in I^S.3\) . except for the case I — 2M and aj = (— l)-' a, of symmetric vortices of family 
(ii) in \S.4^ , except for the case L = M, and of zero-parameter asymmetric vortices of family (iv) 
in J^S.h]) . Asymmetric vortices of family (Hi) in \S. ,5)) can not be continued to the domain e G 0(0), 
except for the only solution \S.26]) . 

Hypothesis 3.8 There exists a one-parameter (modulo gauge transformation) continuation to the 
domain e £ 0(0) of the family of asymmetric vortices yd. 2b]) with < < tt, which includes the 
one-parameter continuation of the two exceptions of Proposition \S.1I\ 

We will not be proving Hypothesis 13 . 81 bevond the second-order reductions. Instead, we use Propo- 
sition to study stability of persistent localized modes of the discrete NLS equation (|2.4|l . 



4 Stability of discrete vortices 

The spectral stability of discrete vortices is studied with the standard linearization: 
where A G C and (a„^m,6„.m) e solve the linear eigenvalue problem on {n,m) G Z^: 

(l ^ Ojijryi 4^n.ni^n.m ^ (^n+l,m '^n— l,m '^n.m+l ^n.rn—l) — i^^n.rai 

~<^n,m'^n,™ + ~ 2|0„^mp) &n,m — £ (&n+l,m + ^n-l,m + &n,m+l + &n,m-l) = —i^bn,m- 

The stability problem (|4.2(l can be formulated in the matrix-vector form: 

Hxp = iXaxp, (4.2) 

where (£ ft x ft consists of 2-blocks of (a„_m, &n,m)"^, is defined by the linearization operator 
H2.12|l . and cr consists of 2-by-2 blocks of 

C -•) 

The discrete vortex is called spectrally unstable if there exists A and i/^ G $7 x $7 in the problem 
H4.2|l . such that Re(A) > 0. Otherwise, the discrete vortex is called weakly spectrally stable. By 
using the Taylor series (|2.17|) . the linearized operator H is expanded as follows: 

H = H^"^ +eH^^^ +e^H^^^ +0{e^), (4.3) 

where Ti^^^ is defined in H2.13|l . while the first-order and second-order corrections are given by 

/ ^(0) ,{i) ,{o) 7(1) ,{o) ,{i) \ /in 

~" .(0)^(1) ^(0) -t-«^(0) ^(1) -(d+i,o + d-i,o + do,+i+do,-i) 



and 

, ,a(0) ^(2) ,(0) 7(2) ,(0) ,(2) \ / ,(1) |2 .(1)2 
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It is clear from the explicit form (|2.13|l that the spectrum of H'^^^if = 7y has exactly N negative 
eigenvalues j = ~2, N zero eigenvalues 7 = and infinitely many positive eigenvalues 7 = +1. The 
negative and zero eigenvalues 7 = — 2 and 7 = map to N double zero eigenvalues A = in the 
eigenvalue problem aTi^^^tp = iXxp. The positive eigenvalues 7 = +1 map to the infinitely many 
eigenvalues A = ±z. Since zero eigenvalues of aTi.^^^ are isolated from the rest of the spectrum of 
(tH(°\ their splitting can be studied through regular perturbation theory |30j . On the other hand, 
if the discrete vortex solutions for (n, m) e decays sufficiently fast as \n\ + \m\ — s- 00, the 
continuous spectral bands of aTi. are located on the imaginary axis of A near the points A = ±i, 
similarly to the case (j>n,m = for {n,m) e [3T]. Therefore, the infinite-dimensional part of the 
spectrum does not produce any unstable eigenvalues Re(A) > in the stability problem (|4.2() with 
small e G 0(0). We shall consider how zero eigenvalues of 7^^°^ and aTi^^^ split as e ^ for solutions 
of the nonlinear equations 1)2. 4|l . which are categorized by Propositions 13. II and 13. 71 



4.1 Splitting of zero eigenvalues in the first-order reductions 

The splitting of zero eigenvalues of Ti. is related to the Lyapunov-Schmidt reductions of the nonlinear 
equations (|2.11() . We show that the same matrix M 1 , which gives the Jacobian of the bifurcation 
functions g^-^\6), defines also small eigenvalues of TC that bifurcate from zero eigenvalues of TY^"' 
in the first-order reductions. 

Lemma 4.1 Let the Jacobian matrix Mi have eigenvalues fJ'^j^\ 1 < J < N. The eigenvalue 
problem Tlip = 7y has N eigenvalues in e Q 0(0), such that 

lim^ = /4^\ l<j<N. (4.4) 

Proof. We assume that there exists an analytical solution 0(£) of the nonlinear equations H2.11|l . 
The Taylor series of 0(f) is defined by H2.17|l . By taking the derivative in e, we rewrite the problem 
(|2.11|l in the form: 

np'ipie) + enM^)+nsCl){e) =0, -0(e) = 0'(£), (4.5) 

where the linearization operator (|2.12|l is represented as 7i = Hp + eHs ■ Using the series H2.17|l and 
1)4. 3|l . we have the linear inhomogeneous equation: 

=0. (4.6) 

Let ej{0), j — 1, ...,7V be eigenvectors of the kernel of Ti'^*'^. Each eigenvector ej{0) contains the 
only non-zero block i(e*^J, —e~^^')'^ at the j-th position, which corresponds to the node {n,m) on 
the contour Sm- It is clear that the eigenvectors are orthogonal as follows: 

(e,(0),e,(0)) ==2<5,,,, l<i,3<N. (4.7) 

Let ej{9), j — 1,...,N be generalized eigenvectors, such that each eigenvector ej{6) contains the 
only non-zero block (e*^^ , e~*^J at the j-th position. Direct computations show that 

n^°^ej{e) = 2icrej(6»), 1 < j < iV. (4.8) 
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The limiting solution H2.5|l can be represented as follows: 

N 

By comparing the inhomogeneous equation (|4.t)|l with the definition H2.15|l of the bifurcation function 
g{0) and its Taylor series (|2.16() . we have the correspondence: 
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Consider a perturbation to a fixed point of g'^^'(0*) = in the form 6 = 6^, + ec, where c 
(ci, C2, cn)^ G R^- It is clear that 



N 



= + e^ce, + 0(e2), e,(0) = e, - ec,e, + 0(6^), 

i=l 

where 0*-°^ — (0,), = 6^(0*), and — ej{9^). By expanding the bifurcation hmction 
g^^\0) near = 0^, we define the Jacobian matrix M.i: 

gf\0)=gfUe{M,c)^ + O{,''), 

where 

1 " 1 ^ 

(-^ic), - - ^(e„H.e,)Q- -c,^(e„Hse,). (4.9) 

1=1 1=1 
On the other hand, the regular perturbation series for small eigenvalues of the problem TL^p = ^ip 
are defined as follows: 

= + 6(^(1) + eV'^^ + 0(e3), 7 = + £2^2 + 0(e3), (4.10) 

where ip^'^^ — 'Ylfj^i'^j^v according to the kernel of 'H^^\ The first-order correction term ^p^^'^ 
satisfies the inhomogeneous equation: 

7^(0)^(1) +^(1)^(0) ^^^^(0)^ (4.11) 
Projection to the kernel of Ti'*'-' gives the eigenvalue problem for 71: 



i^(e„7i(i)e.)c, = 7ic,. (4.12) 

i=l 

By direct computations: 

^ N N 

- 2 E (e^'^-eO = cos{e, - 6,+i) + cos(0, - 0,_i) = ^ (e„H«e,) , 

1=1 i=l 

such that the limiting relation follows from and with TY^^) = T-^p^ +'Hs. m 

We apply results of Lemma EITl to the solutions of the first-order reductions, which are described in 
families (i)-(iv) of Propositions l3.1l and l3.7l The numbers of negative, zero and positive eigenvalues 
of 7W1 are denoted as n{Mi), z{Mi) and p{Mi) respectively. These numbers determine the 
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numbers of small negative and positive eigenvalues of Ti for small positive e. They are predicted 
from Lemmas 13.41 13.51 and 13.61 

For family (i), we compute the parameter Ai in Lemma f3. 41 as Ai — (— 1)^"'(2Z — TV), where / is 
defined in Proposition l3.1l By Lemma we have n(A^i) ~ N — l — l, z{A4i) — 1, and = I 

for < / < 2M - 1 and n{Mi) = N ~ I, z{Mi) = 1, and p{Mi) ^ I - 1 ior 2M + 1 < I < AM. 
In the case of super-symmetric solitons for I = 2M, by Lemma [3.61 we have n{A4i) = 2AI — 1, 
z{Mi) = 2, and p{Mi) ^ 2M - 1. 

For family (ii), by Lemma IXBl we have n{Mi) = 0, z{Mi) = 1, and p{Mi) = N— lioTl<L< 
M - 1 and n{Mi) = N - 1, z{Mi) = 1, and p{Mi) = for M + 1 < L < 2M - 1, where L is 
defined in Proposition 13.11 The case of super-symmetric vortices L = 2M can only be studied in 
the second-order reductions, since Aii = 0. 

The family (iii) is represented by the only one-parameter solution 13.26|l for each L = M. By 
Lemma 15^ we have n{Mi) = 2M - 1, z{Mi) = 2, and p{Mi) 2Af - 1 for < 6* < tt, excluding 
the case of super-symmetric vortices — ^ but including the case of super-symmetric solitons 6 = 
and 9 = 7:. 

The family (iv) is characterized by the value of cos 6** ^ 0, L 7^ M, and 1 < n < — 1, 
n ^ 2M, specified in Proposition 13.11 The parameter Ai in Lemma 13.41 is computed as Ai = 
(— l)^^"(cos0,)^^^(2n — N), such that z{Aii) = 1 in all cases. In the case cos0, > 0, by Lemma 
13.41 we have n{Mi) = N — n — 1 and p{Mi) = n for 1 < n < 2M — 1 and n{M.x) = N — n and 
p{Mi) = n — 1 for 2M + l<n<iV — 1. In the opposite case of cos6'» < 0, we have n(A^i) = n 
and p{Mi) = iV — n — 1 for 1 < n < 2M — 1 and n{Mi) = n — 1 and p{Mi) = N — n iov 
2M + l<n<N -I. 

The splitting of zero eigenvalue of Ti is related to splitting of zero eigenvalues of crTi in the stability 
problem (|4.2|) . 

Lemma 4.2 Let the Jacobian matrix Aii have eigenvalues /ij^^ 1 < J < N. The eigenvalue 
problem Tixp = iXaip has N pairs of eigenvalues Xj m e € 0(0), .such that 

lim^==2^^.^\ l<j<N. (4.13) 

e— »o e ■' 

Proof. The regular perturbation series for small eigenvalues of crTi are defined as follows: 

V> = -0(")+7i,/,(i)+e-i/,(2)^g^^(3)^Q(^2)^ A = x/iAi+eA2 + ex/iA3 + 0(e2), (4.14) 

where, due to the relations 14.7|l and (|4.8|l . we have: 

V.(°)=5:c,e„ V'i) = |^c,e„ (4.15) 

according to the kernel and generalized kernel of <tT-1^^\ The second-order correction term ■0^^-' 
satisfies the inhomogeneous equation: 

7^(0)^(2) + 7^(1)^(0) = ix.aip^^) + iX^cTtP^o) . (4.16) 
Projection to the kernel of Ti'-*'-' gives the eigenvalue problem for Ai: 

Mic=^c, (4.17) 
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where c = (ci, C2, cat)-^ and the matrix A^i is the same as in the eigenvalue problem H4.12|l . As 
a result, the relation (|4.13|) is proved. ■ 



The numbers of negative, zero and positive eigenvalues of Aii, denoted as n{A4i), z{A4i) and 
p{A4i), are computed above. Let ri, zi, and ii be the numbers of pairs of real, zero and imaginary 
eigenvalues of the reduced eigenvalue problem (|4.17|l . We consider stability of families (i)-(iv), 
described in Propositions 13 . II and \3 . 71 

For family (i), we have ii ^ N - I - 1, zi = 1, and n = / for < ? < 2Af - 1; ii = iV - / - 1, 
zi = 2, and n = / - 1 for I = 2M; and ii = iV - zi = 1, and ri = / - 1 for 2M + 1 < / < A^, 
where / is defined in Proposition 13. II 

For family (ii), we have ii — 0, zi = 1, and ri = — 1 for 1 < L < Af — 1; = 0, zi = A^, and 
ri = for L = M; and «i = A^ - 1, zi = 1, and pi = for Af + 1 < i < 2Af - 1, where L is defined 
in Proposition 13. II 

For the only one-parameter solution H3.26(l of family (iii) for each L = M, we have ii = 2M — 1, 
zi = 2, and ri = 2Af - 1 for < 6* < tt and 6* 7^ f . 

For family (iv) with cos 9^ > 0, we have ii=A^ — n— l,zi = l, and ri = n for 1 < n < 2Af — 1 and 
ii = A^ — n, zi = 1, and ri — n — 1 for 2Af + 1 < n < N — 1. In the opposite case of cos 9^, < 0, we 
have ii — n, zi = I, and ri = A^ — n — 1 for 1 < n < 2M — 1 and ii = n — 1, zi = 1, and ri = N — n 
for 2M + 1 < n < N - I. 

There are several features which are not captured in the first-order reductions. For super-symmetric 
solitons of family (i), when I = 2A/ but aj ^ (— l)-'a, the additional zero eigenvalue splits at the 
second-order reductions, which leads to an additional non-zero eigenvalues of the stability problem 
l|4.2|l . For super-symmetric vortices of family (ii), when L = M, the matrix A4i = 0, such that 
non-zero eigenvalues occur only in the second-order reductions. Finally, for symmetric vortices of 
family (ii), multiple real non-zero eigenvalues of the first-order reductions, according to the roots of 
sin^ ^ in the explicit solution H3.13|l . split into the complex domain in the second-order reductions. 
These questions are studied next in the reverse order. 

4.2 Splitting of non-zero eigenvalues in the second-order reductions 

We continue the regular perturbation series H4.14|l to the second-order reductions. By using the 
explicit first-order correction term H3.19|l . we compute the explicit solution: 



where the vectors ej±i are obtained from by shifts of non-zero elements of from the node 
{n,m) S Sm to the adjacent nodes {n,m) e 1?\Sm- The third-order correction term ■0'^-' satisfies 
the inhomogeneous equation: 




N ^ N N 



7^(0)^(3) ^ ^(1)^(1) = + i\2a^^^^ + +iXs<JiP^°\ 



(4.19) 



Projection to the kernel of H^^^ gives the extended eigenvalue problem for Ai and A2: 




(4.20) 
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where the matrix £i is defined by 

(Ak.4f^'^-'')' ;.=^:fj' (4.21) 

[0, b - j| / 1 

subject to the periodic boundary conditions. Let 7j be an eigenvalue of the symmetric matrix Mi 
with the Hnearly independent eigenvector cj. Then, 



Ai = ±v/27,, A2 = ; (\ (4.22) 

Since the matrix Ci is skew-symmetric, the second-order correction term A2 is purely imaginary, 
unless (cj,£iCj) = 0. For discrete solitons of family (i), we have sm{9j+i — 6j) = 0, such that 
£1 = and A2 = 0. 

For symmetric vortices of family (ii) with L ^ M, the matrix Aii has double eigenvalues, according 
to the roots of sin^ ^ in the explicit solution H3.13|l . Using the same discrete Fourier transform as 
in the proof of Lemma 13.51 one can find the values of Ai and A2 in this case. 

Lemma 4.3 Let all coefficients aj = cos{9j+i — 6j) and bj = sin(0j+i — 6j), 1 < j < N be the 



same: aj — a an- 



d bj — b. Eigenvalues of the reduced problem \4-20^ are given explicitly: 

Ai = ±V8asin^, A2 = -2i6sin^, l<n<N. (4.23) 

According to Lemma f4. 31 all double roots of Ai for n 7^ n ^ N split along the imaginary 

axis in A2. When a > 0, the splitting occurs in the transverse directions to the real values of Ai. 
When a < 0, the splitting occurs in the longitudinal directions to the imaginary values of Ai. The 
simple roots at n = ^ and n = N are not affected, since A2 = for n = ^ and n = N. 

For asymmetric vortices of family (iii), it follows from the explicit solution H3.15|l that A2 = for 
all roots, so that the splitting of eigenvalues does not occur in the second-order reductions. 

For asymmetric vortices of family (iv) , the value of A2 can not be computed in general. 

4.3 Splitting of zero eigenvalues in the second-order reductions 

We extend results of the regular perturbation series H4.10|l and (|4.14(l to the case Mi = 0, which 
occurs for super-symmetric vortices of family (ii) with charge L — M. It follows from the problem 
H4.11|l with Ml ^ that 71 = and the first-order correction term tp*^^) has the explicit form: 

^ N N 

where the vectors ej±i are the same as in the formula ()4.18|l . The second-order correction term 
t^^^^ satisfies the inhomogeneous equation: 

7^(0)^(2) + 7^(1)^(1) + n^^)^io) ^ ^2<^(o). (4.25) 

Projection to the kernel of TY^"' gives the eigenvalue problem for 72: 

1 1 ^ 

I (e„H«<^(^)) + (e„H(2)e.) =7^0,. (4.26) 
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By direct computations in the three separate cases, one can show that the matrix on the left- 
hand-side of the eigenvalue problem (|4.26|l is nothing but the matrix A^2, which is the Jacobian 
of the nonlinear function g^-^\6), computed for the super-symmetric vortex of family (ii). Let the 
number of negative, zero and positive eigenvalues of M2 be denoted as n{M2), z{M2) and p{M2) 
respectively. For super-symmetric vortices of family (ii), we have n{A42) = 0, z{Ai2) = 2 and 
p{M2) = 2 for M = 1; n{M2) = 1, z{M2) = 2 and p{M2) = 5 for M = 2; and n{M2) = 4, 
z{M2) = 2 and ^(TWa) = 6 for M ^ 3. 

Splitting of zero eigenvalues of a7i is studied with the regular perturbation series H4.14(l . When 
M.I = 0, it follows from the problem (|4.16|l that Ai — 0, such that the regular perturbation series 
H4.14|l can be re-ordered as follows: 

i/; = V'"' + eV'^^' + eV^^ + 0(e3), A = eAi + e^Az + 0(e3), (4.27) 

where 

^(°)=^c,e„ ^«=<^« + |5:c,e„ (4.28) 

and ip^^^ is given by H4.24|l . The second-order correction term i/)*-^^ is found from the inhomogeneous 
equation: 

7^(0)^(2) + 7^(1)^(1) + 7^(2)^(0) ^ ,Aia^(i) + a2C7t/;(") . (4.29) 
Projection to the kernel of TiS^'^ gives the eigenvalue problem for Ai: 

X2C- Ai/:2C + yC, (4.30) 

where c = (ci, C2, cjv)^, the matrix A^2 is the same as in the eigenvalue problem (|4.26|l . and the 
matrix C2 follows from the matrix Ci in the form H4.21|l with sin(0j+i — Oj) = 1, or explicitly: 

« = J - 1, 

{C2\,,^{ -1, 1 = 3 + 1 (4.31) 

subject to the periodic boundary conditions. Since Ai2 is symmetric and £2 is skew-symmetric, 
the eigenvalues of the problem H4.30|l occur in pairs (Ai,— Ai). Computations of eigenvalues of the 
reduced eigenvalue problem (|4.30(l are reported in the three distinct cases: M — 1^ M = 2 and 
M > 3. 

Case M — 1: The reduced eigenvalue problem H4.30|l takes the form of the difference equation with 
constant coefficients: 

-Cj+2 + 2cj - Cj_2 = AjCj + 2Ai (cj+i - Cj_i) , j = 1,2,3,4, 

subject to the periodic boundary conditions. By the discrete Fourier transform, the difference 
equation reduces to the characteristic equation: 




TTn \ 2 

Ai + 2isinyj =0, ?i = 1,2,3,4. 



The reduced eigenvalue problem (|4.30|) has two eigenvalues of algebraic multiplicity two at Ai = —2i 
and Ai = 2i and zero eigenvalue of algebraic multiplicity four. 
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Case M = 2: The reduced eigenvalue problem H4.3()|l takes the form of the system of difference 
equations with constant coefficients: 

-Xj+i + 2xj - x-i^i = Xlx-i + 2Ai {yj - y-j-i) , j = 1, 2, 3, 4, 
Vj+i - 2j/j+2 + Vj-i = Xfyj + 2Ai {x^+i ~ x^) , j = 1, 2, 3, 4, 

where Xj = (^ij-i and yj — C2j subject to the periodic boundary conditions. The characteristic 
equation for the linear system takes the explicit form: 

At-2A^l-(-l)"-8sin2— j + 8 sin^ — (^1 - (-1)" - 2 sin^ — j =0, n= 1,2,3,4. 

The reduced eigenvalue problem H4.30|l has three eigenvalues of algebraic multiplicity four at 
Ai = — V2j, Ai = \/2i, and Ai = 0, and four simple eigenvalues at Ai = \/80 + 8 and 
Ai = iVVSO- 8. 

Case M > 3: Eigenvalues of the reduced eigenvalue problem H4.3U|I with M = 3 are computed 
numerically by using Mathematica. The results are as follows: 

Ai,2 = ±3.68497i, A3,4 = X5.6 = ±3.20804i, Ay.s = ±2.25068i, A940 = A1142 = ±i, 
Ai3,i4 = Aisae = ±0.53991, Ai7,i8,i9,20 = ±0.634263 ± 0.282851i, A2i,22,23,24 = 0. 

Using these computations of eigenvalues, we summarize that the second-order reduced eigenvalue 
problem (|4.3U|) has no unstable eigenvalues A when L = AI ^ 1; a, simple unstable (positive) 
eigenvalue when L = M = 2; two unstable real and two unstable complex eigenvalues when L = 
M = 3. We note that destabilization of the super-symmetric vortex with M = 2 occurs due to the 
center node (2,2), which couples the four even- numbered nodes of the contour ^2 in the second- 
order reductions. Due to this coupling, there exists a simple negative eigenvalue of the Jacobian 
matrix M2 and a simple positive eigenvalue in the reduced eigenvalue problem (|4.3U|I . Similarly, 
destabilization of the super-symmetric vortex with M — S occurs due to the coupling of eight nodes 
of the contour 6*3 with four interior corner points (2,2), (2,M), (M, M), and (M, 2), which result 
in the four negative eigenvalues of the matrix M2 and four unstable eigenvalues in the reduced 
eigenvalue problem H4.30|l . 

We note that if the matrix AI2 would be defined by the formula (|3.27() for any M > 1 (i.e. if all 
nodes inside the contour Sm would be removed by drilling a hole), the eigenvalues of A^2 would 
be all positive and the eigenvalues of the reduced problem H4.30|l would be all purely imaginary, 
similarly to the case M = 1. 

4.4 Additional splitting in the second-order reduction 

For super-symmetric solitons of family (i), when I = 2M but Oj 7^ (—1)%, the Jacobian matrix 
Ml has two zero eigenvalues with eigenvectors po and pi, but the Jacobian matrix M.i-\- eM2 has 
only one zero eigenvalue with eigenvector pp. Therefore, the splitting of zero eigenvalue occurs in 
the second-order reduction. Extending the regular perturbation series H4.1U|I to the next order, we 
find that 71 = for c = pi, and 

(pi,>l2Pl) 
(Pl,Pl) 
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Extending the regular perturbation series l|4.14|l to the second order, we have find that = for 
c = pi, £i = 0, and 

^2 o (Pl.-^2Pl) „ 

A2 = I— — = 272. 

(Pi,Pi) 

Thus, the sphtting of zero eigenvalue in the second-order reduction is the same as the splitting 
of zero eigenvalues in the first-order reductions. A positive eigenvalue 72 results in a pair of real 
eigenvalues A2 , while a negative eigenvalue 72 results in a pair of purely imaginary eigenvalues A2 . 

All individual results on stability of localized modes on the discrete contour Sm are summarized as 
follows. 

Proposition 4.4 Consider the stability problem i4-^ *^ domain e E O{0), associated to the 
families of discrete solitons and vortices of Provositions and \3.7\ The following solutions are 
spectrally stable in the domain e € 0(0).' discrete solitons of family (i) with I = 0,' symmetric 
vortices of family (ii) with the charge M + I < L < 2M — 1; and symmetric vortex of family 
(ii) with the charge L = M — I. All other solutions have at least one unstable eigenvalue with 
Re(A) > 0. 

We note that stability of discrete solitons of family (i) with / = is equivalent to stability of discrete 
solitons in the one-dimensional NLS lattice, which is proved in the first paper [J. When I = 0, the 
limiting solution 12.5|l consists of alternating up and down pulses along the contour Sm, similar to 
Theorem 3.6 in 

We also note that the purely imaginary eigenvalue of A have negative Krein signature, such that the 
Hamiltonian Hopf bifurcations may occur for larger values of e beyond the neighborhood 0(0). The 
number of eigenvalues of negative Krein signature is related to the closure relation for the number 
of negative eigenvalues of the linearized Hamiltonian H (see [T]). We can see from computations 
of families (i)-(iv) that the closure relation is satisfied in all cases, except for the super-symmetric 
vortices of family (ii) with L — M. Indeed, in the case L = M = 1, the Hamiltonian H has four 
negative eigenvalues for e e O{0), while it has two pairs of purely imaginary eigenvalues of negative 
Krein signature and two pairs of zero eigenvalues, which exceeds the allowed negative index of Ti.. 
Derivation of a modified closure relation for vortex solutions of the discrete NLS equations is beyond 
the scope of this manuscript. 



5 Numerical Results 

We perform direct numerical simulations of the discrete NLS equation l|2.1|l in order to examine 
stability of the discrete vortices in the simplest cases M = 1 and M = 2. The results are shown on 
Figures [1161 fsee also summary in Table 1). 

In computations of solutions of the problems H2.4|l and (|4.2(l . we will use an equivalent renormal- 
ization of the problem with parameter: 



This renormalization is equivalent to keeping the diagonal term — 4e0„^„j in the right-hand-side of 
the difference equations 1)2. 4|l . As a result, the discrete solitons and vortices exist typically in the 
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semi-infinite domain e > 0. The anti-continuous limit is not affected by the renormahzation since 
e « e for small e. 

On the figures the top left panel shows the profile of the vortex solution for a specific value 
of £ by means of contour plots of the real (top left), imaginary (top right) modulus (bottom left) 
and phase (bottom right) two-dimensional profiles. The top right panel shows the complex plane 
A = \r + i\i for the linear eigenvalue problem H4.2(l for the same value of e. The bottom panel shows 
the dependence of small eigenvalues as a function of e, obtained via continuation methods from the 
anti-continuum limit of e = 0. The solid lines represent numerical results, while the dashed lines 
show results of the first-order and second-order reductions. 

Figure ^ show results for the super-symmetric vortex of charge i = 1 on the contour Sm with 
M — 1. In the second-order reduction, the stability spectrum of the vortex solution has a zero 
eigenvalue of algebraic multiplicity 4 and two pairs of imaginary eigenvalues A « ±2ei. These pairs 
split along the imaginary axis beyond the second-order reductions and undertake Hamiltonian 
Hopf bifurcations for larger values of e upon collision with the continuous spectrum (only the first 
bifurcation at e « 0.38 is shown on Fig. It follows from Fig. ^that the zero eigenvalue splits 
along the imaginary axis for larger values of e, beyond the second-order reductions. 

Figure [3 show results for the symmetric vortex of charge L = 1 on the contour Sm with M = 2. 
There are three double and one simple real unstable eigenvalues in the first-order reductions, but 
all double eigenvalues split into the complex plane in the second-order reductions. The asymptotic 
result 14.23(1 for eigenvalues A — -^eAi -I- £\2 with iV = 8, a = cos(7r/4) and h = sin(7r/4) are shown 
on Fig. 121 in perfect agreement with numerical results. 

FigureOlshows results for the super-symmetric vortex with L = M ~ 2. The second-order reductions 
have a pair of simple real eigenvalues A « ±eVv80— 8, a pair of simple imaginary eigenvalues 
zero eigenvalue of algebraic multiplicity 4 and a pair of imaginary eigenvalues of 
algebraic multiplicity 4 at A ~ ±ieV2. The bottom right panel of Fig. shows splitting of multiple 
imaginary eigenvalues beyond the second-order reductions and also subsequent Hamiltonian-Hopf 
bifurcations for larger values of e. 

Figure 21 shows results for the symmetric vortex with L — 3 and M = 2. The first-order reductions 
predict three pairs of double imaginary eigenvalues, a pair of simple imaginary eigenvalues and 
a double zero eigenvalue. The double eigenvalues split in the second-order reductions along the 
imaginary axis, given by H4.23|l with N ~ 8, a = cos(37r/4) and b — sin(37r/4). The seven pairs of 
imaginary eigenvalues lead to a cascade of seven Hamiltonian Hopf bifurcations for larger values 
of £ due to their collisions with the continuous spectrum. The first Hamiltonian-Hopf bifurcation 
when the symmetric vortex becomes unstable occurs for e w 0.096. 

Zero-parameter asymmetric vortices of family (iv) on the contour Sm with AI = 2 are shown in 
Fig. [SI for L = 1 and in Fig. Elfor L = 3. In the case of Fig. |31 all the phase differences between 
adjacent sites in the contour are 7r/6, except for the last one which is 57r/6, completing a phase 
trip of 2tt for a vortex of topological charge L = 1. Eigenvalues of the matrix A^i in the first-order 
reductions can be computed numerically as follows: 71 = —1.154, 72 = 0, 73 = 0.507, 74 = 0.784, 
75 — 1.732, 76 = 2.252, 77 = 2.957 and 73 = 3.314. As a result, the corresponding eigenvalues 
A ~ ±\/2r/e yield one pair of imaginary eigenvalues and six pairs of real eigenvalues, in agreement 
with our numerical results. The bottom panel of Fig. [S] shows that two pairs of real eigenvalues 
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collide for e « 0.047 and e « 0.057 and lead to two quartets of eigenvalues. 

In the case of Fig. El all the phase differences in the contour are 57r/6, except for the last one 
which is 7r/6, resulting in a vortex of topological charge L = 3. Eigenvalues of the matrix A^i are 
found numerically as follows: 71 = -3.314, 72 = -2.957, 73 = -2.252, 74 = -1.732, 75 = -0.784, 
7g = —0.507, 77 = 0, and 73 — 1.154. Consequently, this solution has six pairs of imaginary 
eigenvalues and one pair of real eigenvalues. The first Hamiltonian-Hopf bifurcation occurs for 
e « 0.086. 

We note that numerical results for asymmetric vortices of family (iii) are not shown. Both analytical 
and numerical analysis of such solutions are delicate problems which are open for future studies. 
However, all such solutions (if they persist) are unstable in the first-order reductions. 

6 Conclusions 

In this series of two papers, we have developed the mathematical analysis of discrete soliton and 
vortex solutions of the discrete NLS equations close to the anti-continuum limit. These solutions are 
relevant to recent experimental and numerical studies in the context of nonlinear optics, photonic 
crystal lattices, soft condensed-matter physics, and Bose-Einstein condensates. 

In the present paper, we have examined the persistence of discrete vortices starting from the anti- 
continuum limit of uncoupled oscillators and continuing towards a finite coupling constant e. We 
have found persistent families of such solutions that include symmetric and asymmetric vortices. 
We have ruled out other non-persistent solutions by means of the Lyapunov-Schmidt method. We 
have subsequently categorized the persisting solutions to discrete solitons, symmetric vortices and 
one-parameter and zero-parameter asymmetric vortices. For persistent solutions, we have derived 
the leading-order asymptotic approximations for small unstable and neutrally stable eigenvalues of 
the stability problem, up to the first-order and second-order corrections. 

We have applied the results to particular computations of discrete vortices of topological charge 
L = 1, i = 2, and L = 3 on the discrete square contours Sm with M — 1 and M — 2. Besides 
particular computations collected in Table 1 and Figs. these results offer a road map on stability 
predictions for larger contours, as well as predictions on how to stabilize the discrete vortices. For 
example, super-symmetric vortices of charge L = M > 2 can be stabilized by excluding the inner 
nodes inside the discrete contour Sm- 

There are remaining open problems for future analytical work. First, it would be nice to extend this 
analysis to three-dimensional structures, such as the discrete solitons, vortices, or vortex "cubes" 
(see 123]). Second, other types of contours can be studied in the two-dimensional lattice, such as 
the diagonal square contours which would include the "vortex cross" or the octagon (see |22|)- Per- 
sistence of one-parameter asymmetric vortices including the super-symmetric vortices has not yet 
been proved in the present paper and it leaves space for future work. Finally, the closure relation for 
negative indices of linearized Hamiltonian associated with the discrete vortex solutions must be de- 
rived and applied independently of the current studies. While conceptually, the methodologies and 
techniques presented here can be adapted to the problems mentioned above, actual computations 
of the higher-order Lyapunov-Schmidt reductions become technically involved. 
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Figure 1: The (super-symmetric) vortex cell with L = M = 1. The top left panel shows the profile 
of the solution for e = 0.6. The subplots show the real (top left), imaginary (top right), modulus 
(bottom left) and phase (bottom right) fields. The top right panel shows the spectral plane (Ar, Xi) 
of the linear eigenvalue problem (|4.2(l . The bottom panel shows the the small eigenvalues versus e 
(the top subplot shows the imaginary part, while the bottom shows the real part). The solid lines 
show the numerical results, while the dashed lines show the results of the second-order reductions. 
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Figure 3: The super-symmetric vortex with L = M ^ 2 for e = 0.1. The bottom right panel is an 
extension of the bottom left panel to larger values of e. Remarkable agreement of the theoretical 
predictions (dashed lines) with the numerical results (solid lines) can be observed for small values 
of £. 
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Figure 6: Same as Figure 1 but for the asymmetric vortex with L — 3 and M — 2 for e — 0.1. 
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